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ABSTRACT 

The term 'hypermachine' denotes any data processing device (theoretical or that can be implemented) capable of 
carrying out tasks that cannot be performed by a Turing machine. We present a possible quantum algorithm for 
If^ , a classically non-computable decision problem, Hilbert's tenth problem; more specifically, we present a possible 

^^ ' hypercomputation model based on quantum computation. Our algorithm is inspired by the one proposed by Tien 

^^ , D. Kieu, but we have selected the infinite square well instead of the (one-dimensional) simple harmonic oscillator 

^^ ' as the underlying physical system. Our model exploits the quantum adiabatic process and the characteristics of 

(Ih . the representation of the dynamical Lie algebra su(l, 1) associated to the infinite square well. 

.^ ■ 

Keywords: Hypercomputation, computability, adiabatic quantum computation, infinite square well, dynamical 

[~^ I Lie algebra su(l, 1), Hilbert's tenth problem 



1. INTRODUCTION 



(N 
> 

f^ ' Once upon on time, back in the golden age of 

, -^ . the recursive function theory, computability was an absolute. 

^— ^ ' Richard Sylvan and Jack Copeland.^ 

^^ ' The hypercomputers, according to Copeland and Proudfoot,'^ compute functions or numbers, or more generally, 

(-H , solve problems or carry out tasks, that cannot be computed or solved by a Turing machine (TM) . The hypercom- 

^ i' putation theory rejects the idea of an absolute computability (i.e. Turing computability), detached from logical, 

mathematical, physical or biological theories. Notwithstanding, the up coming of an academic community that 
^ ' works the concept of hypercomputation* , and notwithstanding the proliferation of theoretical hypercomputation 

S ' models,'^ the possibility of real construction of a hypermachine is controversial and is still under analysis. 

O^' At first glance, we could think that the possibility of a hypercomputation model would be a refutation of 

the widely accepted Church- Turing thesis, which is usually interpreted as the identification of the naturally 
calculable functions with the TM-computable functions; but actually the existence of hypercomputation models 
KS ' refute the thesis M, which identifies the functions calculable by a machine with the TM-computable functions.^ 

S . In other words, the existence of the current proposed hypercomputation models coexists with acceptance of the 

Church- Turing thesis. 

The existence of supported proposals of hypercomputation in quantum mechanics is ample, '^ nevertheless this 
is not the case for those based on quantum computation. Considering the fact that quantum computation extends 
beyond its "standard" model (i.e. quantum Turing machines or quantum circuits) and it includes proposals such 
as continuous,^ adiabatic,^ and/or holonomic quantum computation,^ among others; this article presents the 
construction of a hypercomputation model from quantum computation. Our model is based on the one proposed 
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by Tien D. Kieu,^"^" but we have selected the infinite square well (ISW) instead of the (one-dimensional) simple 
harmonic oscillator (SHO) as the underlying physical system. Our model exploits the quantum adiabatic process, 
and due to this change, our model is supported by some characteristics of the dynamical Lie algebra su(l, 1), 
associated to the ISW; instead of the dynamical Lie algebra Weyl-Heisenberg gw^j, associated to the SHO. 

2. KIEU'S ALGORITHM 

A Diophantine equation is of the following form 

D{xi,...,xk)^0, (1) 

where D is a polynomial with integer coefficients. In present terminology, Hilbert's tenth problem may be 
paraphrased as: Given a Diophantine equation of type (1), we should build a procedure to determine whether or 
not this equation has a solution in non- negative integers N. From the concluding results obtained by Matiyasevich, 
Davis, Robinson, and Putnam, we know that, in the general case, this problem is algorithmically insolvable or 
more precisely, it is TM incomputable.^^ 

The hypercomputability of Kieu's algorithm is due to the fact that this algorithm solves Hilbert's tenth prob- 
lem. From its preliminary version, ^^ Kieu has presented different refinements to his algorithm^"^"' i3-i6t_ i^ g^^_ 
dition, with basis in his algorithm, Kieu has reformulated Hilbert's tenth problem in strictly mathematical terms, 
indicating a possible way of solution through the mathematical analysis and the theory of infinite-dimensional 
operators. ^^ 

Kieu's algorithm is constructed beginning with the SHO characteristics. For the SHO with Hamiltonian 

\ (2) 

= a^a + 1/2 , 

the occupation- number states | n) , the action of the creation a^ and annihilation a operators on occupation- 
number states, their commutation relations [a, a^, [a, a], and [a^ , a^, the occupation- number operator N and the 
coherent states | a) are given by 

{|n)|nGN} , (3) 

a|0)=0, a|n) = V^|n-l) , a^ | n) = \^H~T | n + 1) , (4) 

[a, a'] — 1 , [a, a] = [a' , a'] ~ , (5) 

N ^a'fa , (6) 

2 °° n 

l«)-e"^E4il") ' "^C , (7) 

n—O ^ 

With the basis for the SHO, Kieu's algorithm® is shown in Table (1). 
Kieu has indicated the following characteristics of his algorithm®"^°: 

1. The algorithm is probabilistic such as are quantum algorithms in general. 



2. The infinite-dimension Hamiltonian (2) acts on some Fock space whose orthonormal basis is (3), and 
satisfies that | n) ~ ^"'J- \ 0) . 



3. The codification of (1) is made by the occupation-number operator (6) which is diagonal on the basis (3), 
and whose eigenvalues are the non-negative integers. 



^In the Discussion List FOM - Foundations of Matliematics (www.cs.nyu.edu/pipermail/fom/) some cliaracteristics 
of this algorithm have been analyzed and discussed. 



Table 1. Kieu's hypercomputational quantum algorithm. 



Given a Diophantine equation with k unknowns of type (1), Kieu provides the following quantum algo- 
rithm to decide whether this equation has any non- negative integer solution or not: 

1. Construct a physical process in which a system initially starts with a direct product of k coherent 
states I '0(0)), and in which the system is subject to a time-dependent Hamiltonian H^it) over the 
time interval [0, T], for some time T, with the initial Hamiltonian Hi and the final Hamiltonian H^, 
given by 

k 

IV'(0))=(g)|a.) , (8) 

i=l 

HA{t) = (1 - t/T)Hi + (t/T)H^ , (9) 

k 

Hi = '^(al - a*) {ai- ai) , (10) 

Hu ^ iD{Ni, . . . , Nk)f . (11) 

2. Measure through the time-dependent Schrodingcr equation 

idt\^{t))=HA{t)\^{t)), forte [o,r] 

the maximum probability to find the system in a particular occupation-number state at the chosen 
time T, 

Pmax{T) = max |(V'(r) I rii,. ..,nfe)| , 

= \WT)\ {n})f, 

where | {n})^ (which is a direct product of k particular occupation-number state, {^j^x | 'n-i)o) provides 
that maximal probability among all other direct products of k occupation-number states. 

3. If Pmiix{T) < 1/2, increase T and repeat all the steps above. 

4. If 

Pmax(T) > 1/2 (12) 

then I {n})^ is the ground state of Hd (assuming no degeneracy) and we can terminate the algorithm 
and deduce a conclusion from the fact that 



ffD I {"•})o = iff the equation (1) has a non-negative integer solution. 



4. The application of the adiabatic theorem to obtain from the ground state (8) associated with the eigenvalue 
zero of (10), the ground state | {n})o of (11). 

5. The finiteness of the run time T to obtain with a high probability the ground state | {n}{T))Q of Hu. 
Although the value of T is not calculable a priori, its finiteness is guaranteed because the ground state of 
Ha{sT) for < s < 1 is non-degenerate, and this ground state never crosses with any other state during 
the adiabatic regime. 

6. The halting criterion (12) is established by the maximum peak of density of probability associated with 



the initial state (8), where for any n and a 



|2 _ „|a|Ha| 



2n 



n] 
< 1/2 . 

In addition, it is necessary to verify that the probabihty of any excited state cannot be greater than 1/2 
at anytime. 

3. HYPERCOMPUTATIONAL QUANTUM ALGORITHM A LA KIEU 

Although the physical referent used by Kieu for the construction of his algorithm was the SHO, parting from the 
description on Table (1), it can be observed that its constitutive elements, come forth from the dynamical algebra 
Weyl-Heisenberg gvvfl associated with the SHO. The dynamical algebra gw^ satisfies the commutation relations 
(5). Beginning with its generating elements a,a^ and 1, it is possible to factor the Hamiltonian from the SHO 
as pointed out by (2) and it is possible to construct the occupation- number operator (6). The coherent states 
(7) correspond to the eigenstates of annihilation operator a, in other words, a\a) = a|a). Furthermore, the 
algebra gw^j has an infinite-dimensional irreducible representation wherein the action of its generators on basis 
(3) is given by (4). This representation makes posible that the eigenvalues of operator (6) are the non- negative 
integers and its eigenstates are the basis (3), in other words N \n) = n\n). 

To carry out the construction of our algorithm a la Kieu, we have selected as physical referent the ISW, which 
has an associated dynamical algebra different to the one used by Kieu. For a particle with mass m trapped inside 
the infinite well < x < ttI, the Hamiltonian operator Tjisw ^^^j ^j-^g energy levels Ej^"^ are^^ 

fc2 j2 fc2 fc2 

rrlSW _ -2 "• " "" fISW _ "- „f„ , o\ (Tq\ 



where the action of H ° on basis (3) is 



//ISW I . ^ rfSW I 



Due to the spectral structure of the ISW, the dynamical algebra associated with it, is the Lie algebra su(l, 1).^^ 
This is a three-dimensional algebra that satisfies the commutation relations 

[K^ , K+] = Ks , [K^ , K3] = 2X_ , [K+ , X3] = -2K+ , 

where operators K^, K^ and K^ are called creation, annihilation and Cartan operators, respectively. The algebra 
su(l, 1) admits an infinite-dimensional irreducible representation where actions of K^,K^ and Kj, on basis (3) 
are 



K+\n) = v/(n + l)(n + 3)|n+l) , if_ | 0) = , X3 I ") - (2n + 3) | n> 



(14) 



K-\n) = ^Jn(n + 2)\n-l) , 
With basis in the algebra su(l, 1), the Hamiltonian (13) is rewritten as 

and a new number operator N^^^ is given by 

iVisw^(l/2)(i^3-3) , (16) 

where A^^^w | ^^^ _ 'fi\n), in other words, the operator A^^^w pj-^gcnts the same spectral characteristics as the 
occupation-number operator (6). 



Due to the dynamical algebra associated, the Barut-Girardello coherent states | z) , z e C, for the ISW are 
eigenstates of annihilation operator K- ^^ 

|z)= '^' y '^ --\n) , (17) 

where Iv{x) is the modified Bessel function of the first kind. 

With the presented elements, Kieu's algorithm is rewritten in the following way. Instead of replacing each 
one of the variables of (1) by (6) to construct (11), these can be replaced by (16) to obtain 

H]^^={D{Nl^^,...,Nl''^))' . (18) 

Due to (18), it is necessary to construct a new initial Hamiltonian Hj^^ from creation and annihilation (14) 

operators of su(l, 1) 

fc 

ijfw = 5](if+.-z;)(if_,-z.) , (19) 

which has associated with the eigenvalue zero, the ground state I V'^^^(O)), constructed from the coherent states 

(17) 

fe 

|V'''^(0))=(g)|z.) . (20) 

Finally, from (18) and (19) the Hamiltonian (9) takes the form 

Hf^it) = (l - I) Hf-^ + (i) iflfw . (21) 

Some observations with respect to the new algorithm: 

1. The infinite-dimension Hamiltonian (15) acts on some Fock space too, whose orthonormal basis is (3), and 
satisfies that I n) — , ^ ^' I 0). 

' ' ^«!(«+2)!/2 ' ' 

2. The Hamiltonians (9) and (21) arc unbounded operators, therefore it is necessary to use a version of 
the adiabatic theorem for unbounded operators,^" for both algorithms. However, this fact is not very 
significant because the algorithm operates more in the infrared part of the spectrum than in the ultraviolet 
part. 

3. The conditions required for the finiteness of the run time T are satisfied by (21) adapting Kieu's arguments 
for algebra gw-fl* to algebra su(l, 1). 

4. In order to satisfy the halting criterion (12), it is necessary to select Zi values such that \zi\ > 1.6 for the 
construction of (20), since according to the density of probability associated to the states (17) for any n 
and \z\ > 1.6 

\{z\n)\'- 1^1' 



2/2(2 kl) 
<l/2 , 

In the simulations of our algorithm, ^^ this argument has been numerically confirmed. 



5. Unlike the algebra gw-H, the codification of (1) in the algebra 5u(l, 1) could be directly carried out with 
the diagonal opeator K-^. Apparently, due to (14) this codification change, changes the problem to be 
resolved, given that instead of establishing if (1) has, or doesn't have non-negative integers solutions, the 
only thing that could be established is if (1) has or doesn't have solutions of the form 2n + 3, with n e N. 
Nevertheless, this problem is equivalent to Hilbert's tenth problem parting from the construction of an 
finite system of Diophantinc equations. ^^ 

6. Unlike the algebra gv\m, the algebra su(l, 1) admits different types of coherent states (Barut-Girardello, 
Klauder, Perelomov, etc.).^^ The advantage of a type of coherent state over another one would be in 
computational complexity issues, rather than in computational power issues. 

7. The recent equivalence between the adiabatic computation and the "standard" quantum computation^^ 
does not generate any contradiction with the hypercomputational characteristics of our algorithm (or Kieu's 
algorithm), given that our algorithm is one of adiabatic quantum computation over infinite-dimensional 
spaces, and the demonstration of equivalence indicated is for the finite-dimensional case. 

8. The Lie algebra su(l, 1) is the dynamical algebra associated with different physical referents such as the 
infinite cylindrical wells, the Poschl- Teller potentials, quantum optics systems with SU(1, 1) symmetries, 
among others. Therefore, in principle, it is possible to select one of these referents as an underlying physical 
system of our hypercomputational quantum algorithm. ^^ 

4. NON- (TURING MACHINE) HYPERCOMPUTATION 

Within the context of hypercomputation, a hypcrmachine also capable of simulating a universal Turing machine 
is called a super-TM, otherwise it is called a non-TM.^^ The possible universality of our hypercomputation 
model based on the ISW would be established by its capability to generate a set of quantum gates, such that 
any unitary transformation C/(2"), that is, any quantum gate that operates upon n-qubits, can be approximated 
with sufficient exactness by a quantum circuit that is only made of a finite number of gates of this set. 

For example-^-, for the controUed-NOT gate {CNOT'^), the transformation carried out upon the 2-qubit 
canonical basis {| 00) , | 01) , | 10) , | 11) }, is given by 

|x,,)^^^^^^|x,x©,)^ . (22) 

For the ISW with the Hamiltonian independent from the time (13), the states of the system evolve according 
to the Schrodinger equation solution of stationary states^^ 



|n°°(t)) = t/°°(t)|n°°(0)) 



.^ISWj 



|n-(0)) , 
where U°°{t) is the unitary evolution operator and | n°°(0)) = ^^ c„ \n)°° . The U°°{t) matrix elements are 



(23) 



C/„°^(i)=exp(-z^^^|^)^„, , (24) 



where 5np is the Kronecker delta. 



From different choices of the t parameter in (24), and different qubits coding, it is possible to find the 
evolutions corresponding to different quantum gates. In order to build the C'NOT°° gate, based on the normalized 



* Henceforth, we will use the convention of a superindex over the operators (quantum gates) and over the states (qubits) . 
This superindex will denote the dimension of Hilbert space upon which the operators act or upon which the states are 
defined. 



eigenvectors of (22), we code the basis for a 2-qubit i \ 00) , | 01) , - — '-y= — —, - — '—7= — — > in the canonical basis 
(3) by 

\ r,r,\i codification i „i oo lr,i\4 codification i ^i oo 

1 00) > 1 0) , I 01) > I 2) , 

I 10)^1 11)" codification |4)°° + |1)°° |11)"-|10)" codification |4)°°-|1)°° (^5) 



V2 V2 ' \/2 V2 

^ we obtain that CNOT°° ^ U°° (^^ 



According to (24), in t = ^^^ we obtain that CNOT^ = U°° ( 2Hi^ ^ ^^ere the CNOT°° matrix elements 



Q^oc CATOT^lQ^oo ^ I 2^oc CiVOT-^ I 2^c 

\00 IjxOO |-|\00 |i\00 |-i\00 



are {CNOT°°)np — (— l)"(5„p, and the transformation carried out upon the coded basis (25) is given by 
|0)°°- 

|4)°° + |1)°° CNOT^^ |4)°°-|1)°° |4)°°-|1)°° CNOT^^ |4)°° + |1) 

^/2 ' V2 ' V2 ' \/2 

Given that the universal minimal sets of quantum gates contain gates of 1-qubit and the CNOT gate, it is 
not possible to carry out a codification of the gates of 1-qubit over the same basis used for the CNOT"^ gate, 
therefore, it is not possible to obtain the universality of our hypercomputation model parting from compound 
universal sets via gates that act upon a different number of qubits^. A possibility would be to obtain the 
TM-universality parting from unitary sets of universal gates, such as the Toffoli or Fredkin gates. This will be 
investigated elsewhere. 

5. CONCLUSIONS 

It is quite surprising that a quantum computation model over such a simple physical system like the infinite square 
well has hypercomputational characteristics as it was shown. The success obtained by choosing an underlying 
physical system different from the one selected by Kieu opens the possibility of obtaining new hypercomputation 
models supported by the quantum computation. 

Kieu's algorithm does not depend on any intrinsic peculiarity of the simple harmonic oscillator nor on its 
dynamical algebra gw-fl- It is possible to consider that in principle any finite-dimensional dynamical algebra 
that admits an infinite-dimensional irreducible representation, and that it admits the formation of attainable 
coherent states in a known quantum system, is a good candidate to establish an algorithm d la Kieu. 
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